QCD factorization for B PP 



Dongsheng Du^'^, Deshan Yang^ and Guohuai Zhu^ ^ 
^ CCAST (World Laboratory), P.O.Box 8730, Beijing 100080, China 
^ Institute of High Energy Physics, Chinese Academy of Sciences, P.O.Box 918(4), Beijing 

100039, China ^ 
(May 18, 2001) 

Abstract 

In this work, we give a detailed discussion for QCD factorization in- 
volved in the complete chirally enhanced power corrections in the heavy 
quark limit for B decays to two light pseudoscalar mesons, and present 
some detailed calculations of radiative corrections at the order of a^. 
We point out that the infrared finiteness of the vertex corrections in the 
chirally enhanced power corrections requires twist-3 light-cone distribu- 
tion amplitudes (LCDAs) of the light pseudoscalar symmetric. However, 
even in the symmetric condition, there is also a logarithmic divergence 
from the endpoints of the twist-3 LCDAs in the hard spectator scat- 
tering. We point out that the decay amplitudes of -B ^ PP predicted 
by QCD factorization are really free of the renormalization scale depen- 
dence, at least at the order of a^. At last, we briefly compare the QCD 
factorization with the generalized factorization and PQCD method. 
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I. INTRODUCTION 



The study of B decays plays an important role in understanding the origin of CP viola- 
tion and physics of heavy flavor. We expect that the parameters of the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix in the standard model, for instance, the three angles a, (3 and 7 in 
the unitary triangle, can be well-determined from B decays, especially from the charmless 
non-leptonic two-body B decays. Experimentally, many B experiment projects have been 
running (CLEO, BaBar, Belle etc.), or will run in forthcoming years (BTeV, CERN LHCb, 
DESY HeraB etc.). With the accumulation of the data, the theorists will be urged to gain 
a deeper sight into B decays, and to reduce the theoretical errors in determining the CKM 
parameters from the experimental data. 

In the theoretical frame, the standard approach to deal with such decays is based on 
the low-energy effective Hamiltonian which is obtained by the Wilson operator product 
expansion method (OPE). In this effective Hamiltonian, the short- distance contributions 
from the scale above /i ~ have been absorbed into the Wilson coefficients with the 
perturbative theory and renormalization group method. The Wilson coefficients have been 
evaluated to next-to-leading order. Then the main task in studying non-leptonic two-body 
B decays is to calculate the hadronic matrix elements of the effective operators. However, we 
do not have a reliable approach to evaluate them from the first principles of QCD dynamics 
up to now. 

Generally, we must resort to the factorization assumption to calculate the hadronic ma- 
trix elements for non-leptonic B decays, in which the hadronic matrix element of the ef- 
fective operator (in general, which is in the form of current- current four-quark operator) 
can be approximated as a product of two single current hadronic matrix elements; then it 
is parameterized into meson decay constant and meson-meson transition form factor. The 
most popular factorization model is the Bauer-Stech-Wirbel (BSW) model [jT|. In many 
cases, BSW model achieves great success, which can predict the branching ratios of many 
modes of non-leptonic B decays in correct order of magnitude. This factorization assump- 
tion does hold in the limit that the soft interactions in the initial and final states can be 
ignored. It seems that the argument of color-transparency can give reasonable support to 
the above limit. Because h quark is heavy, the quarks from h quark decay move so fast 
that a pair of quarks in a small color-singlet object decouple from the soft interactions. But 
the shortcomings of this simple model are obvious. First, the renormalization scheme and 
scale dependence in the hadronic matrix elements of the effective operators are apparently 
missed. Then the full decay amplitude predicted by BSW model remains dependent on the 
renormalization scheme and scale, which are mainly from Wilson coefficients. In past years, 
many researchers improved the simple factorization scheme and made many remarkable 
progresses, such as scheme and scale independent effective Wilson coefficients effective 
color number which is introduced to compensate the 'non-factorizable' contributions, etc. 
Furthermore, some progresses in nonperturbative methods, such as lattice QCD, QCD sum 
rule etc. allow us to compute many non-perturbative parameters in B decays, such as 

the meson decay constants and meson-meson transition form factors. Every improvement 
allows us to have a closer look at the B nonleptonic decays. 

Except for the factorization approximation, another important approach has been ap- 
plied to study many B exclusive hadronic decay channels, such as i? ^ Dvr, tttt, t^K etc. 
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This is PQCD method In this method, people assumes that B exclusive hadronic 

decay is dominated by hard gluon exchange. It is analogous to the framework of perturba- 
tive factorization for exclusive processes in QCD at large momentum transfer, such as the 



calculation of the electromagnetic form factor of the pion ||10[. The decay amplitude for 
B decay can be written as a convolution of a hard-scattering kernel with light-cone wave 
functions of the participating mesons. Furthermore, in Ref. the Sudakov suppression 
has been taken into account. 

Two year ago, Beneke, Buchalla, Neubert, and Sachrajda (BENS) gave a QCD factoriza- 
tion formula in the heavy quark limit for the decays B ^ [TT]]. They pointed out that the 
radiative corrections from hard gluon exchange can be calculated by use of the perturbative 
QCD method if one neglects the power contributions of AgcD/^b- This factorization for- 
mula can be justified in case that the ejected meson from the b quark decay is a light meson 
or an onium, no matter whether the other recoiling meson which absorbs the spectator quark 
in B meson is light or heavy. But for the case that the ejected meson is in an extremely 
asymmetric configuration, such as D meson, this factorization formula does not hold. The 
contributions from the hard scattering with the spectator quark in B meson are also involved 
in their formula. This kind of contribution cannot be contained in the naive factorization. 
But it appears in the order of a^. So they said that the naive factorization can be recovered 
if one neglects the radiative corrections and power AQco/^b suppressed contributions in the 
QCD factorization, and the 'non-factorizable' contributions in the naive factorization can 
be calculated perturbatively, then we do not need a phenomenological parameter A^^-^-^ to 
compensate the 'non-factorizable' effects any more [p!2|-[T^. 

This QCD factorization (BBNS approach) has been applied to study many B meson 
decay modes, such as -B — >■ D^*^tt~ |T5| , p!6| , tttt, ttK UTTHTsI and other interesting channels 
pil|^3[. Some theoretical generalizations of BBNS approach have also been made, such as 
the chirally enhanced power corrections |lT8| , p!9| , |2^ , ^ from the twist-3 light-cone distribution 
amplitudes of the light pseudoscalar mesons. In this work, we will take a closer look at this 
issue. This work is organized as follows: Sect. II is devoted to a sketch of the low energy 
effective Hamiltonian; in Sect. Ill, we will give a detailed overview of QCD factorization, 
in which some elaborate calculations are shown, especially for the chirally enhanced power 
corrections; Sect. IV is for some detailed discussions and comparison of BBNS approach to 
the generalized factorization and PQCD method; we conclude in Sect.V with a summary. 



II. EFFECTIVE HAMILTONIAN FIRST STEP FACTORIZATION 



B decays involve three characteristic scales which are strongly ordered: mw ^ ^fe 3> 
Aqcd- How to separate or factorize these three scales is the most essential question in B 
hadronic decays. 

With the operator product expansion method (OPE), the relevant \AB\ = 1 effective 
Hamiltonian is given by ||26|| : 
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lq=u,c \ k=3 
- Vt {C7^{fl)Q7^{n) + C8G(Ai)(58G(/^))] + h.C, 



(1) 
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where Vq = VqbV*^(foY b ^ d transition) or Vq = V^feVg^ (for b ^ s transition) and Ct(/i) are 
the Wilson coefficients which have been evaluated to next-to-leading order approximation 
with the perturbative theory and renormalization group method. 

In the Eq.(l), the four-quark operators Qi are given by 

Q" = {Uaba)v^A{qf3U^)v^A Q\ = {Caba)v-A{qi3Ci3)v-A 

Q2 = {uJ)p)v-A{q(3Uo)v-A Ql = {Cabf3)v-A{qi3Ca)v-A 

Q3 = iqaba)v-AEiqf3qfs)v-A Q4 = iqt3ba)v-AEiqaq't3)v-A 



Qb = iqaba)v-AY.iqi3qi3)v+A Qg = iqf3ba)v-AY.iqaqf3)v+A 

q' q' 

Q7 = U^aba)v-AEeqiq,3q'p)v+A Qs = U^(3ba)v-AE eq'{q^q'/^)v+A 



(2) 



<?' q' 



Q9 = \{.qMv-~AT.eq^{.qpqp)v~A Qw = \{qpba)v-Alleq'{q^qp\ 
q' q' 



V-A 



and 



Q,^ = —m,q^a^^-{l + ^^)b^F^,, Qg^ = ^ + 75)^/3^^., {q = dois). (3) 
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with Q\ and Q\ being the tree operators, — Qe the QCD penguin operators, Q-j — Qio 
the electroweak penguin operators, and Qy-y, Q^g the magnetic-penguin operators. 

In this effective Hamiltonian for B decays, the contributions from large virtual momenta 
of the loop corrections from scale /i = 0{mij) to mw are attributed to the Wilson coefficients, 
and the low energy contributions are fully incorporated into the matrix elements of the 



operators [^. So the derivation of the effective Hamiltonian can be called "the ffist step 
factorization" . 

To evaluate the Wilson coefficients, we must extract them at a large renormalization scale 
(for example = 0{my/) in the standard model) by matching the amplitude of the effective 
Hamiltonian {A(.ff) to that of the full theory {Afuii)-, then evolve them by the renormalization 
group equations from the scale /i = 0{mw) to the scale /i = 0{mi,). It should be noted 
that the extraction of the Wilson coefficients by matching does not depend on the choice 
of the external states, if we regularize the infrared (and mass) singularities properly ||26|| . 
All dependence on the choice of external states only appears in the matrix elements {Qi), 
and is not contained in Cj. So Ci only contains the short-distance contributions from the 
region where the perturbative theory can be applied. But for the matrix elements (Qi), the 
long-distance contributions appear, and are pro cess- dependent. 

Several years ago, the perturbative corrections to the Wilson coefficients in SM have been 



evaluated to next-to- leading order with renormalization group method |2^. As we know 



the Wilson coefficients are generally renormalization scheme and scale dependent. So, in 
order to cancel such dependence, we must calculate the hadronic matrix elements of the 
effective operators to the corresponding perturbative order with the same renormalization 
scheme and at the same scale, then we can obtain a complete decay amplitude which is free 
from those unphysical dependences. 

III. QCD FACTORIZATION FOR B ^ PP 

After "the ffist step factorization", the decay amplitude for B hih2 can be written as 
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A{B ^ hih^) = ^t;,Q(/x)(/ii/i2|Qi(/i)|B) , (4) 

i 

in which, as mentioned in the previous section, the contributions from the large scale mw 
down to mf, has been separated into the Wilson coefficients Ci{fi). The remaining task is to 
calculate the hadronic matrix elements of the effective operators. But for the complexity of 
QCD dynamics, it is difficult to calculate these matrix elements reliably from ffist principles. 
The most popular approximation is factorization hypothesis, in which the matrix element of 
the current- current operator is approximated to a product of two matrix elements of single 
current operator: 

(hihmB) ~ (/i2|J2|0)(/ii|Ji|B) . (5) 

Obviously, under this approximation, the original hadronic matrix element {Qi{fi)) misses 
the dependence of the renormalization scheme and scale which should be used to cancel the 
corresponding dependence in the Wilson coefficients Ci{fi). A plausible solution to recover 
this scale and scheme dependence of (Qi) is to calculate the radiative corrections. In one- loop 
level, they can be written as |1^,|T^J2^: 



(Q) = [i + ^m, + ^rhe] ■ (Q)TVee. (6) 

Here ihg and liie represent the one loop corrections of QCD and QED respectively. Then 
one takes 

(/il/i2|Q^|B)Ttee ^ (/i2 | ^2 1 0) (/^l | Jl | B) . (7) 

Therefore, the scheme and scale dependence of (Qi) which are expressed in the form of riis 
and rile is recovered. But in quark level, liis and liie usually contain infrared divergences if 



we take the external quarks on-shell |]28|. To remove or regularize the infrared divergence, 
the conventional treatment is to assume that external quarks are off-shell by — p^. But this 
introduction of the infrared cutoff — results in a gauge dependence of one-loop corrections. 
So how to factorize the infrared part of the matrix elements is a very subtle question. But 
maybe this question would get an important simplification in the case that the final states 
of B meson decay are two light mesons. 

Two years ago, Beneke, Buchalla, Neubert and Sachrajda proposed a promising QCD 
factorization method for B — >■ vrvr. They pointed out that in the heavy quark limit 3> 
AqcDi the hadronic matrix elements for B can be written in the form 

(vrvrigiB) = {'K\.h\Q)Wi\B) ■ [1 + + ^(AqcdM)]. (8) 

Obviously, the above formula reduces to the naive factorization if we neglect the power 
corrections in KQcn/fT^b and the radiative corrections in a^. They find that the radiative 
corrections, which are dominated by hard gluon exchange, can be calculated systematically 
with the perturbative theory in the limit rrtb — > cxd, in terms of the convolution of the hard 
scattering kernel and the light-cone distribution amplitudes of the mesons. This is also 
similar to the framework of perturbative factorization for exclusive processes in QCD at 
large momentum transfer, such as the calculation of the electromagnetic form factor of the 



pion |jTO|. Then a factorization formula for B ^ tttt can be written as ||11|| : 
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W>{q)mB{v)) = F^'^^iq') C dxTl{x)^^{x) + didxdyT'\E,,x,y)^B{i)^.{x)^M- 

JO Jo 

(9) 

We call this factorization formalism as QCD factorization or the BBNS approach. In the 
above formula, ^it{x) are the leading-twist wave functions of B and pion mesons 

respectively, and the T/'^^ denote hard-scattering kernels which are calculable in perturbative 
theory. At the order of as, the hard kernels T^'^^ can be depicted by Fig.l. Figures 1(a)- 1(d) 
represent vertex corrections. Figs 1(e) and 1(f) penguin corrections, and Figs 1(g) and 1(h) 
hard spectator scattering. 

In the heavy quark limit, both pions are energetic. The pion ejected from b quark decay 
moves so fast that it can be described by its leading-twist light-cone distribution amplitude. 
The qq pair in the ejected pion is produced as a small-size color dipole. Consequently, the 
ejected pion decouples from the soft gluons at leading order of KQCD/i^h- Of course, only the 
cancellation of soft gluons is not enough to make the factorization hold, it is necessary that 
the qq pair also decouples from the collinear gluons. Both the cancellations of soft gluons 
and collinear gluons guarantee that the hard kernel T/ is of infrared finiteness. Contrast to 
the pion ejected from b quark weak decay, the recoiling pion which picks up the spectator in 
B meson can not be described by its leading- twist light-cone distribution amphtude(LCDA), 
because the spectator is transferred to the recoiling pion as a soft quark. Here Beneke et 
al. take the point of view that the form factor F^^'^ cannot be calculated perturbatively. 
If we attempt to calculate the form factor within the perturbative framework, by the naive 
power counting, we find that the leading twist LCDA of pion does not fall fast enough to 
suppress the singularity at the endpoint where the quark from h decay carries almost all 
momentum of the pion. It indicates that the contributions to form factor are dominated by 
the soft gluon exchange |]TB|. This point of view can be justified also from the calculation 



of the form factor F^^^ by using light-cone sum rule (LCSR) , in which the dominated 
contribution to F^^'^ comes from the region where the the spectator quark is transferred as 
a soft quark to the pion. So the transition form factor survives in the factorization formula 
as a nonperturbative parameter. However, when the spectator quark in B meson interacts 
with a hard gluon from the ejected pion, the recoiling pion can be also described by its 
light-cone distribution amplitude. This hard spectator scattering is missed in the naive 
factorization, but calculable in the perturbative QCD at the leading power in Kqcd/tti}). 
So with this factorization formula, the remaining hard part of the hadronic matrix element 
{Qi) from the scale about has been factorized into the hard scattering kernel, and the 
long distance contributions are absorbed into the transition form factors and the light-cone 
wave functions of the participating mesons. Thus this is the "final factorization" for the 
two-body nonleptonic charmless B decays. 

An explicitly technical demonstration of the above argument has been presented in one- 
loop level in Refs. p^JT6[] . For B Dtc, this QCD factorization has been proved to two-loop 



order |T6|. In the literature, the ejected pion is represented by its leading twist light-cone 
distribution amplitude(LCDA). However, since the mass of b quark is not asymptotically 
large, in particular, some power corrections might be enhanced by certain factors, such as 
the scale of chiral symmetry breaking /i^ = m^/(m„ + m^) ~ 1.5 GeV, and have significant 
effects in studying B two-body nonleptonic charmless decays. So, in this manner, the chirally 
enhanced power corrections must be taken into account. Accordingly, describing the ejected 
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pion by its leading twist LCDA is not enough, the two-particle twist-3 LCDAs must be 
taken into account. Below, we will show the elaborate results of QCD factorization in these 
two cases. For illustration, we take tt+tt" as an example, but the result is easily 

generalized to the cases that the final states are the other light pseudoscalars. 



A. Leading-twist Distribution Amplitude Insertion 

When inserting leading-twist LCDA of the light pseudoscalar, in the heavy quark limit, 
the quark constituents of the ejected pion can be treated as a pair of coUinear massless 
quark and antiquark with the momentum uq and uq respectively (g is the momentum of 
the ejected pion and we take g as a hard light-cone momentum in calculation, u = 1 — u), 
because that the contributions from the transverse momenta of the quarks in ejected pion 
are power suppressed [|l^ . 



1. Vertex Corrections 

Now we move on to the explicit one- loop calculation of the diagram Figs. l(a)-l(d) for 
B TTTT. For illustration, we write down the one-gluon exchange contribution to the B^ —>■ 
TT+TT^ matrix element of the operator Q2 = {uabf3)v-A{df^Ua)v-A = (<^a&o)F-A(w/3'U/3)v-^- 



d^k 1 



(27r)4 - ky[{p - kf - ml] 

x(7r+|M,7^(l - 75)^757° W -^)7,.(1 - 75)(?/-^ + m6)7,6,|5°), (10) 

{Q2){b) =9s-r^ / du(f){u 



4 N Jo ^-^v-^y {2n)^ k^{uq - ky[{p - ky - ml] 
X (7r+|n,7^(l - 75)(n^ - ^)7"!^757^(1 'l^W-H m,)7«6.|^d), (H) 

X (7r+|M,7«(2/ - ^ -^)7^(1 - 75)^757aW + ^)7m(1 - 75)&.|S°), (12) 



A N Jo ' J {2TiYk'^{uq + kf{p-q-kf 

x(7r+|n,7,(2/- ^ -^)7^(1 - 75) (^^ + ^) 7a ^757^(1 - 75)&.|^°), (13) 

When we calculate the vertex corrections in the leading power of KqcD/^hi not only 
ultraviolet divergence emerges but infrared divergence does also. Infrared divergence arises 
from two regions where the virtuality of the loop k is soft or collinear to the momentum of the 



pions. In Ref. [|16[, the authors gave an explicit cancellation of soft and collinear divergence 
in vertex corrections for B Dt^ in eikonal approximation. Figures 1(a), 1(b) and 1(c), 1(d) 
cancel the soft divergence; 1(a), 1(c) and 1(b), 1(d) cancel the collinear divergence. For 
B TTTT, the cancellation is similar except that the collinear divergence also arises from 
the region where k is collinear to the momentum of the recoiling pion. So Figs. 1(c), 1(d) 
cancel not only part of soft divergence but also part of collinear divergence. Below, we give 
an explicit calculation of the Feynman diagrams Figs. l(a)-l(d) to show the cancellation 
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of the infrared divergences. In order to regularize the infrared divergence, there are two 
choices for us. One is the dimensional regularization (DR) scheme, in which the infrared 
divergence can be regularized into the pole terms l/((i — 4). In contrast to the dimensional 
regularization of ultraviolet divergence, the infrared divergence arises when d < 4, instead 
of d > 4 in the case of the ultraviolet divergence. So the dimension d in regularization for 
infrared divergence must be set to be greater than 4. This is a subtle point, but it will 
not cause any ambiguity in our calculation because the infrared part and ultraviolet part 
can be safely separated. The other method to regularize the infrared divergence is the well- 
known massive gluon (MG) scheme, in which the infrared divergence is handled by replacing 
l/k"^ by — m'^) in the gluon propagator. Similar scheme has been applied in earlier 

computation of the radiative corrections for — > e'Ve^^i-, in which the massless photon is 
replaced by a massive photon. In addition, in our latter calculation, there are also several 
schemes in treating 75, the most popular two are the naive dimensional regularization (NDR) 
scheme and the 't Hooft-Veltman renormalization (HV) scheme. Both have been applied 
to calculate the Wilson coefficients p6[. In this work, if there is no specification, the NDR 
scheme is always applied in our calculations for its simplicity ^. 

After a straightforward calculation in DR scheme and using the corresponding Feynman 
parameter integrals listed in Appendix C, we obtain 



m^a) = ^%(vr-|47'^(l-75K|0)(vr+|^^7,(l-75)&/.|5°) 



4:71 N 

X / du(l){u) 



r(l — a) I rrih 
(47r)« y /i 



1 , , M u 

- - 7b + In 477 + 2 In — + 1 + 

e rUb I — u 

2a 



Inu 



1 2(lnM-l) 



+ 



+ ln^M - 2Li2(l - 



4 In M + 5 + 



21nM 



u 



u 



47r 

X / du(j){u) I —4 



1 , , , u 11 -u 
- - 7ij + In 47r + 2 In — + — + 

t VTib 4 \ — u 



Xtvu 



+ 



r(l — a) f nib 
(47r)" \ fi 



2a 



1 2(lnM - 1) , 2 . / Ix 

— + ^ + In^ M - 2Li2(l --)-41nM + 6 + 



21nM 



u 



u 



7£; + In 47r + 2 In — H In(-u) 

e rrib 4 



X J du(f){u) <! -4 



+ 



r(l — a) /nib 
(47r)" I yu 



2a 



2 2(ln(-n) - 2) 



+ 



TT 



+ 10 - 41n(-M) + In^(-M) 

3 



Such choice does cause a scheme dependence in the matrix elements. However, when we choose 
the Wilson coefficients in the same scheme as for the matrix elements, the final full decay amplitude 
is free of scheme-dependence. 
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ducpiu) 



- - 7s + ln47r + 21n — + 1 - \\ii-u) 



rfi 



(4vr)» 



2a 



2 , 2(ln( 



+ 10 



TT 



- 41n(-M) + In^ 



-u] 



In above, we have set = 4 + 2a (a > 0) in regularizing the infrared divergence. Then, 
after summing over all four diagrams, we find that all pole terms in 1/a are really cancelled 
before we integrate over the momentum fraction variable u. So after modified minimal 
subtraction {MS), we get 



{Qfl {a)+{b)+{c)+{d) 

^^(vr-|47^(l - 75)«„|0)(vr+|^,7.(l - l.)h,m 
1 

duipiu 

nib 







-1^ 



121n-^ + 



u 



1 



Inw — 



4m 



u 



1 



u 



In u + 4 ln(— m) — ln(— m) 



- (in^n-ln^w) +2 (^Li2(l 
/ 2 In M 2 In M 



1, 



u 



-Li2(l 



1, 

u' 



u 



u , 



+ {\n^{-u) - In^ 



-u 



Assuming that the light-cone distribution amplitude 
equation can be simplified as follows: 



(18) 



{u) is symmetric, then the above 



{Qt) {a)+{b)+{c)+{,d) 

(7r-|J„7'^(l - 75)n„|0)(7r+bi^7^(l - i,)hp\Bl) 



X 



as Cf 

4^1V 
1 

du(f){u 





-1^ 



121nA + 3LJ!^lnw 



rrib 



1-u 



(19) 



It is easy to check that the above equations are consistent with the results in previous 
works. Actually, with the MG scheme, we get the same results as that by using the DR 
scheme. 

With Eqs.(|l^,|l^), we can compute the vertex corrections no matter the LCDA 0(m) is 
symmetric or asymmetric. This is very important in principle. For instance, when kaon is 
ejected from b quark decay, we must take the contributions from the asymmetric part of 
LCDA of kaon into account, although the contributions from the asymmetric part are very 
small numerically |[T9|. 



2. Penguin Corrections 

There are two kinds of penguin corrections. One is the four quark operators insertion 
[Fig. 1(e)]; the other is the magnetic penguin insertion [Fig. 1(f)]. The first kind is 
generally called BSS mechanism. In generalized factorization, BSS mechanism plays a very 
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important role in CP violation because it is the unique source of strong phases. But in 
generalized factorization, the virtuality of the gluon or photon is ambiguous; usually one 
varies k"^ around This variation does not have an important effect on the branching 

ratios, but it does for CP asymmetries. In QCD factorization, this ambiguity is rendered 
by taking the virtuality of the gluon as k'^ = {p — uqY = uml/2. When treating penguin 
contractions, one should be careful that Fig. 1(e) contains two kinds of topology, which is 
depicted in Fig. 2. They are equivalent in 4 dimensions according to Fierz rearrangement. 
However, since penguin corrections contain ultraviolet divergence, we must do calculations 
in d dimensions where these two kinds of topology are not equivalent Below we give 
an explicit calculation of Q4 or Qq penguin insertions for tt+tt^ which belong to the 

second topology. Fig. 2(b): 



//O \twist— 2 
\V4,6/(e) 



X 



^^s/^^'^ ^ du 0(M)(7r+|Mi7„^757^(l - 75)&i|^d) 
^ f d^k -Tr[(//-^-m,)7"(^ + m,)7'^(lT75)] 



a. Cf 



[(/-A;)2-m2][P-m2]/2 



=p—uq 



9afi 



X 



E 



1/1 , . A 

;:( 1E + ln47r) 

6 e 



dt t{l - t) In 

Jo 



m? - t{l - t)P - it 



l=p~uq 



(20) 



After MS subtraction and using the equations of motions, we get the finite result 

-twist— 2 



(Q4,6) 



X 



E 



a.. C 
4 



. — ^(tt-I ^,7^1 - 75)^dO)(vr+|n,7,(l - 75)&,|5°) 



lnA_4 



du (f){u) dt t{l - t) ln(s„ - t(l - t)u - ie] 
Jo Jo 



(21) 



where Sg = rrig/ml. The first topology. Fig. 2(a), for example, Ql penguin insertion for 
— > 7r+7r~, is similar to the results of the second topology, Fig.2(b), except that there is 
an extra factor —2/3: 

as Cf 



{QT) 



■twist— 2 



A-K N 
2 4 



{'K-\da^{l - 75)w.|0)(7r+|M,7^(l - 75)fe,|5°) 



+ - In — - 4 

3 3 rrif, 



/•I 

du (t){u) I dt t{l — t) ln(sc — t(l — t)u — ie) 
Jo 



For the magnetic penguin insertion, it is the easiest calculation of the radiative correc- 
tions. The result of Qsg insertion for B^ 7i~^ti~ is 



//O \ twist— 2 



as_C 
'An N 
as Cp 
AttW Jo 



'^^fTrrrib I du (j){u)^{'K^\ui'j"^'j5apak'^{l + 'y5)bi\B'^) 
1 2(f){u 



k=p—uq 



du 



U 



-(7r-|rf,7^(l - 75)|0)(vr+|^i,7,.(l " 75)1^^") 



(23) 
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3. Hard scattering with the spectator 



Hard spectator scattering [Fig. 1(g) and (h)] is completely missing in the naive factoriza- 
tion. But in QCD factorization, it can be calculated in the perturbative QCD, and expressed 
by a convolution of the hard kernel T^^ and the LCDAs of mesons. At the leading power of 

AqcD/nT'h, both of the hght pseudoscalars from the B meson decay can be represented by 
their leading twist LCDAs. So after a straightforward calculation, we obtain this contribu- 
from the operator Q2 insertion. 



tion for 



/^•u\ twist— 2 
\V2/(g) + (ft) 



64 m 



X 



9s 



du dv (f)B{0 'Pi'^) 'Pi'^) 



+ mB)757"^i757p(l - 75)^2757a/d7p(l - Ts)] 



272 



Tr[(j/ + mB)757"^i757p(l - 75)/«7a^2757p(l - 75)] 



+ 



• , r2CF ['.Mi) 



r du^ r 

Jo u Jo 



k=^p—vqi 

id=m2-k 

k=^p—vqi 
l„=k~uq2 



dv 



(24) 



B. Chirally enhanced corrections — Twist-3 LCDAs insertion 

It has been observed that QCD factorization is demonstrated only in the strict heavy 
quark hmit. This means that any generahzation of QCD factorization to include or partly 
include power corrections of AQco/nT-b should rcdcmonstratc that factorization still holds. 
There are a variety of sources which may contribute to power corrections in l/m^; examples 
are higher twist distribution amplitudes, transverse momenta of quarks in the light meson, 
annihilation diagrams, etc. Unfortunately, there is no known systematic way to evaluate 
these power corrections for exclusive decays. Moreover, factorization might break down 
when these power corrections, for instance, transverse momenta effects, are considered. This 
indicates that one might have to give up such an ambitious plan that all power corrections 
could be, at least in principle, incorporated into QCD factorization order by order. One 
might argue that power corrections in B decays are numerically unimportant because these 
corrections are expanded in order of a small number Aqcn/'mb — 1/15. But this is not 
true. For instance, the contributions of operator to decay amplitudes would formally 
vanish in the strict heavy quark limit. However it is numerically very important in penguin- 
dominated B rare decays, such as interesting channels B — > nK, etc. This is because 
Qq is always multiplied by a formally power suppressed but chirally enhanced factor r^. — 
2mp/m5(mi+m2) ~ C>{^), where mi and m2 are current quark masses. So power suppression 
might probably fail at least in this case. Therefore phenomenological applicability of QCD 
factorization in B rare decays requires at least a consistent inclusion of chirally enhanced 
corrections. 

Chirally enhanced corrections arise from the two particle twist-3 light-cone distribution 
amplitudes, generally called 0p(x) and (pai^)- So when chirally enhanced corrections are con- 
cerned, the final light mesons should be described by leading twist and twist-3 distribution 
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amplitudes. Then it is crucial to show that factorization really holds when considering twist- 
3 distribution amplitudes. The most difficult part is to demonstrate the infrared finiteness 
of the hard scattering kernels T/. In addition, possible chirally enhanced power corrections 
can also appear in the hard spectator scattering. So, for consistency, we must involve these 
corrections. 



1. Vertex corrections 

When we calculate the chirally enhanced power corrections, contrast to the leading-twist 
light-cone wave function insertion, the cancellation of the infrared divergences in the vertex 
corrections to {V — A){V + A) operator (here it is Q5 or Qj) can not be shown simply by 
the eikonal approximation similar to what has been done at the leading power of h.Qco/'fnb-i 
because the Dirac structure or spin structure of twist-3 light-cone wave functions of the light 
pseudoscalar makes the "on-shell" condition for the external quarks invalid. Thus, to justify 
the cancellation of the infrared divergence in {V — A){y + A) vertex corrections, we must 
give the explicit calculation. As mentioned in the previous subsections, we have two choices 
to regularize the infrared divergence in one-loop calculation. One is the DR scheme; the 
other is MG scheme. Generally, these two schemes are equivalent, for instance, similar to 
what has been done in {V — A){y — A) vertex corrections. However, in the DR scheme, it is 
difficult to extrapolate the twist-3 wave functions of the light pseudoscalar to d dimensions 
properly, although they are well-defined in 4 dimensions. Therefore, we prefer to use the 
MG scheme in our calculation for chirally enhanced corrections to avoid the above possible 
problems. 

In addition, generally we calculate the Feynman diagrams in the momentum space, so 
the correct projection of the light-cone wave functions of the light pseudoscalar in the mo- 
mentum space is necessary. From Appendix B, we find that it is easy to find the proper 
momentum space projection of the leading twist and 0p type twist-3 wave function, but for 
00-, the projection is not very clear. Note that the coordinate in the definition of 4>a by 
the non-local matrix element must be transformed into a partial derivative of a certain mo- 
mentum in the projection of momentum space. However, it is difficult to find the derivative 
which makes the projection only depend on the structure of the light pseudoscalar itself. 
Generally, the momentum which the partial derivative acts on is dependent on the hard ker- 
nel. Therefore, we prefer to compute the Feynman diagrams of the twist-3 wave functions 
insertion, especially (j)„ insertion in the coordinate space. We think that such calculation 
can avoid the ambiguity about how to project the coordinate into the momentum space. 
We recalculate the leading twist insertion by using the same method, and obtain the same 
results as those in the previous sections. Below, we will show how to perform this trick in 
calculation of insertion. For example, let us consider Fig. 3. In coordinate space, we 
have 

' = ^ aT^^ Jdu-^ Jdx,dx,J (^(^(^ 

_fpfipCF 2/", Mu) f d^k d%d*X2 e^("'?-'=+'^)-^2 

~ aT^^ J J (2^ (27r)4 - ml)ll{ll - ml) 
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lb=p-k 



ifpfip Cp 



N 



9s J du 



{2ny {k^ - mml - ml) 



X 



d 



dL 



{7T+\u, 7^(1 + 75)f 7"75^^^.7p(l - 75)(A + m,)7,6,|50) 



iu=k—uq 



lh=p-k 



(25) 



The above trick has been apphed in the calculation of the proper correlation function to 
extract the transition form factor F^^'^ within the frame of the light-cone sum rule 0. So 
within the same frame, we obtain 



(<95)(a) = g: 



/tt/^tt Cp 

4 



du 



d^k 



1 



(2vr)4 (P - m2)/2(/2 _ ml) 



X {0p(n)(7r+|Mi7''(l +75)757"Ma.(1 -75)(A + "^fe)7a&i|5°) 



+ t 



.4>a{u) 



(7r+|M,7''(l + 75)75a;..gV 7' 



(Q5)(6) = ^: 



4 N Jo 



du 



7p(l - 75) (A + mb)^abi\B'^ 

d^k 1 
{2-^)^ {k^ - m^)ll{ll - ml) 



la=uq-k 
lb=p-k 



X 



{(/)p(M)(7r+|Mi7'^(l +75)A7°757^(1 -75)(/fe + mb)7^6i|S°) 



+ « — ^(7r^|ni7^(l + 75) 7 - 



{Qb){c) = g. 



6 



/tt/^tt Cp 



II 



du 



7p(l - 75) (A + rribhabilB'^ 
d^k 1 



lu=k—uq 
lt=p-k 



4 iV Jo "~~J {27i)Hk' - ml)Pjl 
X {0p(u)(7r+|Mj7c,/„7^(l + 75)757a/rf7p(l - 75)&i|i?°) 

0a (m) 



+ ^ 



6 



■(7r+|ui7"A7''(l + 75)75C^p^g''7Q 



{Q5)(d) =g\ 



4 N Jo 



du 



d'^k 



I'd 



7p(l-75)fc.|5°) 



ld=uq-k 



lu=p-q+k 



X {(j)p{u){7T+\ua"^ul^{l + 75)/«7a757p(l " 75)&i|5°) 



+ i 



6 



(vr+|n,77n7''(l+75) W' 



II 



Perform the one-loop integrations: 



lu=k—uq 

lu=p-q+k 



(26) 



(27) 



(28) 



(29) 
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X 



(- -7£; + ln47r + 21n— ) 
+ ^In^ A + ln(-M)lnA - 2lnulnX + - In^ u - Liafl --) + - + ^tt'^ 



+ 



6u 



1 



ln^A + 21n(-'u) In A - 41nuln A + In A 



7 



+ 41n^ii - In^(-w) -21niiln(l-w) + 2Li2(-) + -7r -Inii- 



lnw 



= -2^^(7r-|(i.(l + 75)«^|0)(7r+|«,(l-75)&,|^,°) 



M 6 

■0\ 



u 



X J du l4>p{u) 



-(- -7E + ln47r + 21n— ) 
/o L ' L e mft 

1 11 

+ - W A + Inf-iZ) In A - 2 IniZln A + - W - Liafl ) 

4 2 u 



2 + r 



+ 



6u 

2 . 



^ In^ A + 2 In ( -m) In A - 4 In M In A + In A 



+ 41n^'f2-ln^(-M) - 2 lniZln(l - + 2Li2(4) + -tt^ -IniZ- 

K 6 1 — 



(Q5>(c) = -2^^(7r-|J,(l + 75)^^^|0)(7r+|%(l - 75)&,-|^°) 



1 

X / du 
1 



-(- -7£ + ln47r + 21n— ; 
1 



+ - In' A - [In(-M) - 2] In A + ^ ln"(-M) - In(-M) 



+ 



2 
6m 



3 

2 + y 



In^A- [21n(-M) -3] In A 



TT 



+ 21n(-'u)lnii-ln''u-31n(-w) +3- — 

o 



47r 



X 



du |0p(m 



-(- -7E + ln47r + 21n— ) 



TT 



+ ^ W A - [In(-iZ) - 2] In A + ^ \^{-u) - In(-iZ) + | + „ 



+ 



2 

0a (m) 



In^A- [21n(-iZ) -3] In A 



TT 



+ 21n(-'u)lnM-ln''u-31n(-M) + 3- — 



(30) 



(31) 



(32) 



(33) 



Here A = m^/m^. Prom the above equations, it is observed that, in the case of ^„ distri- 
bution amphtudes, the terms with infrared divergence in vertex correction diagrams cannot 
cancel unless 0o-('^) is a symmetric function: 0o-('") — ^c{^- This is an unexpected result. 
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which means QCD factorization is violated for asymmetric twist-3 hght-cone distribution 
ampUtudcs. This indicates that chirally enhanced corrections can be included consistently 
in the framework of QCD factorization only when twist-3 light-cone distribution amplitudes 
are symmetric. Therefore, in the following, we will implicitly assume a symmetric twist-3 
light-cone distribution amplitude for light pseudoscalar mesons. It is then straightforward 
to show that vertex corrections of {V — A){V + A) operator are completely cancelled after 
summing over four diagrams in the case of 0o- distribution amplitude. The final result of 
{V — A){y + A) vertex corrections, in the condition that the twist-3 LCDA is symmetric, is 

(g5)(a)+(6)+(c)+(d) = 12^^(7r-|(i,(l + 75)«^|0)(7r+|«,(l-75)&,|5°) (34) 



2. Penguin corrections 



In quark level, usually one decomposes the basic QCD vertex —igs'~f'^T°j in penguin 
insertion into the two chiral current couplings —i^gs'j'^Tf^j{l + 75) and —i^gsj^T^j{l — 75); 
then the penguin insertions contribute the same magnitude to the {V — A){V — A) and 
{V — A){V + A) vertex. But in hadron level, this point of view must be examined in 
elaborate calculation. 

For illustration, we give the results of Q4 or Qq penguin corrections to i?^ 7r"'"7r~, 
which belong to the second penguin topology Fig. 2(b), when 4>p{u) is inserted: 

\V4,6;(e) —9s^^ 



X 



E 



d'^k -Tr[(;/ -uj-]^- mg)7"(^ + 771^)7^(1 T Ts)] 
{T.'kY [{p — uq — ky — mV\[k'^ — m^ij) — uqY 



47r iV JO 



X 



-(- 



7E + ln47r)+ C dtt{l-t)\u 
Jo 



11^ 



ml - t{l - t)P - ie 



l=p—uq 



(35) 



After MS subtraction, we obtain 

7r-|(i,(l + 75)«i|0)(7r+|«,(l - 75)&,|5°) 

V fin - 3 du (f>p(u) C dt t(l - t) ln(s„ - t(l - t)u - ie) 
„ I rrih Jo Jo 



OisCp 



'l/l iV 

nl /■! 

-3/ 

q - rrib 

For the first penguin insertion topology. Fig. 2(a), the result is 



X 



(36) 



mTe^ 



— 2— — 
~ 47r 



X 



(7r-| Ji(l + 75)iii|0)(7r+|%(l - 75)&,-|^°) 

/ du <j)p{u) I dt t(l — t)ln(sc — t(l — t)u — ie) 
Jo Jo 



mb 2 



(37) 
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Similarly, when (f)a{u) is inserted, by using the method in the previous subsection, we 
have 



(^4,6) 



4 '^'^^ N Jo 



X 



E 



d f d% Tr[(/,-^ + m,)7"(/, + m,)7''(lT75)] 



o r 
dkn J 



(27r)' 



[{k-kY-mmi-mm 



(38) 



k=p—uq 



After integration and subtraction, 

a., C 



{Q 



4,6/ 



3 mft 



X / du — 
g Jo 6u 

- dt { 2t(l - t) ln(s„ - t(l - - ie) + — ' 

h \ So — t 1 



(39) 



t(l — t)ti — ie ^ 

The magnetic penguin insertion is easier; we write the result of Q%g insertion as follows: 

irt \twist-3 _ OisCp 



X 



£ du |(/-p(^/)^(7r+|^/i7"75(T/jaA:'^(l + 75)^'^|^d> 



+ ^ 



. 4><t{u) 1 



6 4A;2 



{T^^lUinMla] [7^7a](l +75)&i|^°) 



- 7^(^^l«.7"75[^,^ [^,7j(l + 75)&^|^« 



A;2 

c 
)_ 

'47r TV '2 'Jo 



k=p—uq 



^2^^il+ tdu Mf))(^-|J.(i + ^,)^.|o)(7r+|xZ,(l-75)6,-|^S). (40) 



5. ^fard scattering with the spectator 

The chirally enhanced power corrections in hard spectator scattering not only occurs in 
the case of {V — A){V + A) vertex insertion, but also in the case of {V — A){y — A) insertion. 
But in the case of {V — A)(y -\- A) insertion, after a straightforward calculation, we find 
that there will be serious linear divergence at the end points of the LCDAs if the twist- 
3 LCDAs are not symmetric. Because infrared finiteness of the vertex corrections requires 
that the twist-3 LCDAs, especially (pai^), must be symmetric, we shall implicitly assume this 
symmetric condition for the LCDAs in latter computation. So, in this symmetric condition, 
the hard scattering with the spectator vanishes when (V — A){V + A) vertex is inserted. 
However, even in this strict symmetric condition, there is still a logarithmic divergence from 
the end point of the recoiling pion in hard spectator scattering when [V — A){V — A) vertex 
is inserted. For example, 

{Q2){g)+{h) = ^TraJs/^— d^ du — — dv 



(t){v) ^ 2ii^(p^{v) 



V niB Qv"^ 



(41) 
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This means that QCD factorization is broken down. But we can still give a phenomeno- 
logical treatment for this hard spectator scattering. By using the asymptotic form of 0o-('u), 
we find that there is a divergent integral over v: Jq dv (l/v). In Refs. [^,^, the authors 
prefer to introduce a phenomenological parameterization for this logarithmic divergence. 
They take dv (l/v) = \n{mB/AB) + re*^, where r is taken from 3 (realistic) to 6 (conser- 
vative), and the phase 6 from — tt to tt. We shall take similar phenomenological treatment 
in the numerical computation below. 

We notice that the above approach of evaluating hard spectator contribution is naive. 
For instance, the scale of hard spectator contribution should be different from the vertex 
correction contribution. While it seems reasonable to take the scale fi ~ 0{mh) for the ver- 
tex correction diagrams to avoid large logarithm log(/i/mb), a natural choice of the scale 
of hard spectator contribution may be around 0{1 GeV) because the average momentum 
squared of the exchanged gluon is about 1 GeV"^. Another disturbing feature of hard spec- 
tator contribution is that, as pointed out in Refs. [T^J2^, when including the contribution 
of there would appear a divergent integral Jq dv{l/v) even if the symmetric distribution 
amplitude is applied. This divergent integral implies that the dominant contribution comes 
from the end-point region, or, in other words, it is dominated by soft gluon exchange. How- 
ever, the transverse momentum may not be omitted in the end-point region if so, the 
corresponding divergent integral would then be changed to 

f. nv,k 



V J v^ml + 

As an illustration, we do not consider the kx dependence of wave functions (though it is 
certainly not a good approximation); then the above integral is proportional to 

dvdk'^ f dxdy ^^^^ 



v^ml + kj' J X + y 

The above integration converges now; furthermore it is not dominated by end-point con- 
tribution. This illustrates that the treatment of hard spectator diagrams may need further 
discussion. 

There exists "annihilation topology" contributions which may belong to chirally en- 
hanced corrections. In Ref. |2^, the authors have discussed this topic and find that a 
divergent integral [/(rfx/x)]^ will appear. We suspect that this divergence may disappear, 
similar to the hard spectator term, if the effect of transverse momenta can be included. It is 
also possible that "annihilation topology" contributions are really dominated by soft interac- 
tions and thus violate factorization. Due to its complexity, we do not include "annihilation 
topology" contributions in this work. 

C. Final formulas 

With these effective operators, B P1P2 decay amplitudes in QCD factorization can 
be written as 

A{B P,P,) = % E E ^X(^i^2|g.|B)^, (44) 

p=u,ci=l lo 



17 



where Vp is CKM factor, (PiP2|Qj|B) p is the factorized matrix element and is the same as the 
definition of the BSW Largrangian [l|]. Then as an illustration, the explicit expressions of af 
(i = 1 to 10) for 5 — i> vrvr (using symmetric LCDAs of the pion) are obtained. But it is easy 
to generalize these formulas to the case that the final states are other light pseudoscalars. 
Furthermore, we take only part of QED corrections into account in our final formula, in 
particular the QED penguin insertions. Now af for i? ^ tttt in NDR 75 scheme is listed as 
follows ^: 

-~ {c,(ilog ^ + G(s,) - |) + (C, - ^)(| log A + G(0) + G(l) - i) 

+ E {C, + Ce + -e^Cs + -e^Cio){-log^ + G{sg)) + GsCsG\, (48) 

q=u,d,s,c,b 

+ 1^.) log ^ - ^ - + G'isp) + G'^(s,)) 
47r A L 3 m;, 2 3 

+ E (C3 - ^)((l + |a.) log ^ - ^ - + G"(.,) + C^is,)) 

q=d,b ^ ^ TUb 2 6 



+ E (C4 + Ce + -eqCg + -CqCio) ( (1 + -A^) log-^ + G"(s,) + G"^(s,) ) 

q=u,d,s,c,b 2 2 \ 6 rUb J 

+ (^ + A.)C8g}, (50) 
«r = C^r + ^ + ?^%Cs(-F-12), (51) 



A^ 477 A^ 
iV ~ 4^ aT 



^Because of the tedium, we do not calculate the radiative corrections in the HV scheme. However, 
generally, the results in the NDR scheme and HV scheme can be related by a constant matrix 
Afs = fg^uv — rs,NDR [29 1 which is free from the gauge dependence and infrared structure of the 
theory. Thus, in principle, we can obtain the results in the HV scheme just by using Afg. In 
| ]23| , the constant matrix has been applied to obtain the results in the NDR and HV scheme for 
the coefficients Oj which only contain the current-current vertex corrections. But whether we can 
obtain the expression of og or og in HV scheme in a similar way needs further discussion. 
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q=d,b ^ 



q=u,d,s,c,b ^ 



a,-C, + — + ^^C,„F, 



q=u,d,s,c 



Here = 3 is the number of color, Cp = {N'^ — 1)/2A^ is the factor of color, Sg — m 
and we define the other symbols in the above expressions as 

F = -121n-^ - 18 + / + /^^ 

rub 

1 1 
f^j dx g{x)(f){x), Gs^ j dx Gs{x)(f){x), 



1 

G{s) = j dx G{s,x)(p{x), 



1 

G'{s) ^ J dx G'{s,x)(l)p{x), 





where (f){x) [4>p{x), 4>a{x)] is leading twist (twist-3) LCDA of the ejected pion, and 
hard-scattering functions are 



1 — 2x 2 
g{x) — 3— In a; — Siir, Gsix) — - — 

-L X i 



X 



G{s, x) = —4: J du u{l — u) ln(s — u{l — u){l — x) — ie), 



1 

x) — —3 J du u{l — u) ln(s — u{l — u){l — x) — ie). 
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x) = -2 J du u{l - u) ln(s - u{l - - x) - ie) 



u\l-u)\l-x) 



+ du 



s — u{l — u){l — x) — ie 



(63) 



The contributions from the hard spectator scattering [Figs. 1(g), (l)(h)] are reduced to the 
factor f^': 



f 



11 47r2 fjs 



N Ff^^(0)m| 



\i^[dx't^(dy 



X 



l-y Mb 6(1 - yy 



(64) 



There contains a divergent integral in f^^. Here we simply assume that J{dy/y) ~ 
ln{mi,/ Kqcd) (similar to what has been done in Refs. ||T9| , |25[| , though our assumption here 
is certainly an oversimplification). We thus illustrate numerically the scale dependence of 
(tttt) in Table 1. Here we use the asymptotic form of the LCDAs of the light pseudoscalar 
meson which are listed in Appendix A, and the other input parameters are taken as follows 
ig: F^^(O) = 0.33, /b = 0.2 GeV , = 133 MeV, the pole masses rrib = 4.8 GeV, = 
1.4 GeV, the JTS masses mt{mt) = 170 GeV, mb{mb) = 4.4 GeV, m„(2 GeV) = 4.2 MeV, 
md{2 GeV) = 7.6 MeV and Ag^^ = 225 MeV. 



IV. DISCUSSIONS AND GENERAL REMARKS 
A. Color transparency and factorization 

Color transparency gives a clear physics picture of QCD factorization. In the argument 
of the color transparency, a fast-moving small color singlet formed by a pair of qq decouples 
from the surrounding soft gluons. Of course, as mentioned in the previous section, only the 
decoupling with the soft gluons is not enough for a factorization formula, the decoupling from 
the coUinear divergence is also necessary. We really find that both of the requirements can be 
satisfied in the one-loop calculations. So the QCD factorization is guaranteed. Therefore, 
the calculations in the above sections seem to be a one-loop technical manifestation (or 
demonstration) of the color transparency. On the other hand, at the leading power of 
^QCD/mb, the soft or colhnear gluons only "see" the direction of the light meson, but 
are "blind" to the spins of the quark constituents. So the soft or colhnear gluon cannot 
distinguish whether the ejected meson from b quark decay is a light pseudoscalar or a light 
longitudinally polarized vector meson. As a consequence, the cancellation of the infrared 
divergence is universal for B decays to two light mesons, no matter whether the meson 
is a pseudoscalar meson or a vector meson. Therefore, the QCD factorization formula for 
B — > PP at the leading power of Aqcd / f^h is easy to be generalized to B PV and VV . 

Similarly, the color transparency argument can not only be applied to the strong interac- 
tions, but also be generalized to the electromagnetic interactions. When the ejected meson 
is electric neutral, the soft photons also decouple from the fast moving small electric dipole. 
So QED vertex corrections are also of infrared finiteness. But for the case that the ejected 
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meson is charged, QED corrections are infrared divergent, and the infrared divergence must 
be cancelled by the soft photon emission mechanism, which is common in the calculation 
of the radiative corrections for /i~ e~t'^z/e. About this, it is easy to be covered in the 
calculation in the previous section, just replacing the QCD vertex by a QED vertex. This 
can be called a one-loop demonstration for "charge" transparency. 

It should be noted that the above arguments must be based on the condition that the 
ejected meson is in a very compact configuration, then it, as a small color dipole, is disen- 
tangled with the soft gluons. Otherwise, if its size is too large, it is difficult to decouple from 
the soft gluons. For example, the spectator quark in D meson is very soft, and runs around c 
quark like a soft quark cloud, which has a large overlap with the B meson spectator system 
|T6| . As a consequence, the process in which D meson is ejected from b decay is dominated 
by the soft gluon exchange. 



B. The scale dependence 

From the expressions of the QCD coefficients obtained in previous sections, it is found 
that the renormalization scale dependence of the hadronic matrix elements of the effective 
operators is recovered. Apparently, we expect this recovered dependence can cancel the scale 
dependence of the Wilson coefficients Cj. 

With the renormalization group equations for the Wilson coefficients Cj(yu) at leading 



order logarithm approximation |26 



we do find 

^-^a^ = (for i = 1 - 5 and 7, 9, 10) (66) 

dfx 

when we neglect the contributions from higher order of a^. But for ag or ag, some scale 
dependence at the order of still remains. Note that other QCD coefficients (01=1,2,3,4,5,7,9,10) 
are multiplied by the product of the matrix elements of the conserved currents which are 
independent of the renormalization scale; whereas the coefficient og or ag is multiplied by a 
product of the two matrix elements of scalar and pseudoscalar current 

-2(Pi|rf(l + 75)g|0)(P2|g(l-75)b|B). 

which is still of scale dependence. This scale dependence is generally represented by the 
factor 



2ml 



mb(/i)(mi(/i) +m2(/i)) 



after we apply the equations of motion to transform the (5* + P){S — P) matrix elements 
into the type of {V — A){V — A). Here mi and m2 are the current masses of the valence 
quarks in meson Pi. With the renormalization group equations for the running mass of the 
current quark 
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we have 



Consequently, we find 



/^^^(/^) = -6^CFm(/x), (67) 



d ex 

/^^^x(/") = ^'^^CprM ■ (68) 



(<8^x) = (69) 



with the asymptotic form of (paiu) = 6u{l — u). Then, as we expect, the decay amplitude for 
B decays to two light pseudoscalars predicted by QCD factorization is really independent 
of the renormalization scale within the constraint = 1/2. This also can be obviously 
seen from the numerical results of af(7r7r) listed in Table 1. In particular, if we think that 
the results of QCD factorization are reliable and really independent of the renormalization 
scale, maybe Ag- = 1/2 is a strict constraint for the form of 0cr('u). 

It should be noted that the imaginary part of QCD coefficients only arises at the 
order of a^, and depends on the renormalization scale. This dependence would bring some 
uncertainties in determining the CP asymmetries in B decays. Maybe this scale dependence 
of the imaginary part could be cancelled by the results on higher order of as- 



C. Comparison to the generalized factorization and PQCD method 

Comparing QCD factorization approach with the generalized factorization [0,0] and 



PQCD method ||3T[, some interesting comments are in order. 



(i) At the zeroth order of a^, both of QCD factorization (BBNS approach) and generalized 
factorization can reproduce the results of "naive factorization"; at the higher order of as, 
the renormalization scheme and scale dependence for the hadronic matrix elements can be 
recovered from the hard-scattering kernels T/ in BBNS approach and in generalized 
factorization. However, in generalized factorization, uis is from the one-loop calculations 
of quark- level matrix elements. According to Buras et al, quark-level matrix elements are 
accompanied with infrared divergences. To avoid these divergences, one may assume that 
external quarks are off-shell. Unfortunately, it will introduce gauge dependence which is also 
unphysical. But in the BBNS approach, it is different because the external states are all 
physical and can be approximated as on-shell quarks in the leading order of Aqcn/'^b- As a 
consequence, the unphysical gauge dependence does not appear. In the PQCD method, Li et 
al. claim that their method is based on a six-quark system in which the external quarks are 
on-shell, so the gauge invariance of PQCD predictions is guaranteed. The scale dependence 
in PQCD prediction is removed by evolving the Wilson coefficients down to the proper hard 
scale. So no explicit scale dependence is left. 

(ii) In generalized factorization and BBNS approach, the hadronic transition form factors 
are not calculable, and they are dominated by soft gluon exchange and determined only by 
experiments or some non-perturbative approaches such as sum rule, lattice QCD, etc. In 
particular, the above assumption can be also justified in the BBNS approach by naive power 
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counting [|TT] , p!6| . However, in the PQCD method, this naive power counting rule may be 
invahd when the transverse momenta of the quark constituents and Sudakov suppression 
are taken into account. So, Li et al. thought that the hadronic transition form factors can 
be calculated in the PQCD method because h quark is heavy enough and the soft gluon 
exchange is suppressed by Sudakov form factors. This is the essential difference between the 
BBNS approach and the PQCD method. 

(iii) The generalized factorization considers "nonfactorizable" contributions as in- 
tractable. Therefore, one may introduce one or more effective color numbers N^^^ to phe- 
nomenologically represent "nonfactorizable" contributions |T^n]. Furthermore, N^^^ is 
assumed to be universal to maintain predictive power. However, "nonfactorizable" contri- 
bution is really process-dependent. In the BBNS approach and the PQCD method, such 
"nonfactorizable" contributions are indeed calculable in perturbative theory. In consequence, 
N^f^ need not be introduced. 

(iv) As mentioned in the above sections, the strong phases predicted by generalized fac- 
torization are only from the BSS mechanism which is represented by the penguin insertion. 
However, the virtuality of the gluon or photon k'^ in the penguin insertion is ambiguous in 
generalized factorization, and usually it is approximated around m\/2. This brings signifi- 
cant uncertainties for predicting the CP asymmetries for B decays. A particular interesting 
result of the BBNS approach is that strong phases are not only from the BSS mechanism but 
also from the hard scattering, and there are no uncertainties in determining fc^ of penguin 
insertion. However, compared with the real part of the decay amplitude, the imaginary part 
is 0{as) or power KqcD/'mb suppressed and cannot lead to large CP asymmetries, since 
they come solely from hard scattering processes which are only calculable in the heavy quark 
limit. In the PQCD method, there is no such ct^ suppression in the imaginary part of the 
decay amplitude. Thus CP asymmetries predicted by the PQCD method are usually greater 
than the prediction of BBNS approach and generalized factorization. So maybe these dif- 
ferences of prediction for CP asymmetries can be an experimental test for the BBNS and 
PQCD approaches. 

(v) Hard spectator contributions [Fig. l(g)and 1(h)], which are leading power effects in 
QCD factorization, miss out in "naive factorization" and "generalized factorization" . They 
are, however, 0{as) suppressed compared with the leading factorized contributions (the 
hadronic transition form factors). But, in the PQCD method, they are of the same order as 
the form factors. 

(vi) In the PQCD method, penguin contributions receive a dynamical enhancement called 
"Fat Penguin" [§. But in generalized and QCD factorization, they are missed. This en- 
hancement in the PQCD method arises from the strong scale dependence of the penguin 
Wilson coefficients 6*4,65 etc. 

(vii) Final state interactions (FSI) do not appear in the three methods. In QCD factor- 
ization, Beneke et al. point out that the cancellation of the infrared divergences implies that 
the long distance FSI is power suppressed due to the quark-hadron duality. However, this 
point of view is controversial B^, but can be examined by the experimental measurements 



of B ^ KK |3|. 
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D. Limitation of QCD factorization 



QCD factorization formula only holds in the heavy quark limit mt, oo. In the real 
world rUb is only about 4.8 GeV, the validity of the power suppression may be questionable. 
In particular, for several cases, the power suppressed corrections can be numerically large 
||16|| , because the perturbative expansion is in order of as which is not small at the realistic 
scale 0{mh) compared to KqcD/t^b- 

(i) The hard "nonfactorizable" contributions computed by QCD factorization are gen- 
erally small compared to the leading "factorizable" contribution. But when the leading 
"factorizable" contributions are color suppressed, the "nonfactorizable" contribution may 
be larger than the leading results. At the same time, the potentially soft contribution, 
which is formally power suppressed, may be important. For example, in —>■ tt^tt^, any 
perturbative and soft power suppressed contributions can have a significant effect on pre- 
dicting the branching ratio and CP asymmetry. Furthermore, this problem also arises when 
the entire leading power contribution is suppressed by small Wilson coefficients, for exam- 
ple, in B ^ KK; or when the leading power contribution is suppressed by the small CKM 
elements. 

(ii) An important power suppressed contribution is from the higher twist light-cone wave 
functions of the light mesons. The chirally enhanced power correction from the two-particle 
twist-3 wave functions is the most important, and has been partly involved in this work 
except for the annihilation topologies. Other contribution from multiparticle non-valence 
fock state has been proved to be also power suppressed |T^. However, there is no systematic 
way to evaluate it. The author of Ref. proposed a way to evaluate the soft gluon 
exchange contribution from higher twist qqg wave functions within the frame of the light- 
cone sum rule(LCSR). But the accuracy of LCSR is limited due to the quark-hadron duality 
approximation. On the other hand, power correction from transverse momenta needs a 



subtle treatment in the future. In Ref. ||16[, the authors point out that the contribution 
from the transverse momenta might be considered when we evaluate the hadronic matrix 
elements to two-loop order. In this case, it is possible that Sudakov suppression might be 
taken into account as well. 

In summary, up to now, we do not have a systematic way to evaluate many kinds of 
power suppressed corrections for exclusive processes. How to evaluate such corrections in a 
consistent way within the frame of QCD factorization is a potentially interesting work. 



V. SUMMARY 

In this work, we give a detailed discussion for QCD factorization involving the complete 
chirally enhanced power corrections in the heavy quark limit for B decays to two light 
pseudoscalar mesons, and present some elaborate calculations of radiative corrections at the 
order of a^. We point out that the infrared finiteness of the vertex corrections in the chirally 
enhanced power corrections requires twist-3 light-cone distribution amplitudes (LCDAs) of 
the light pseudoscalar symmetric. However, even in the symmetric condition, there is also 
infrared divergence from the end point of the LCDAs in the hard spectator scattering and 
annihilation topology. So the transverse momenta and Sudakov suppression should be taken 
into account. We also point out that the decay amplitude of B ^ PP predicted by QCD 
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factorization is really independent of the renormalization scale, at least at the order of as- At 
last, we briefly compare the QCD factorization to the generalized factorization and PQCD 
method which are generally used in studying B exclusive hadronic decays. 
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Appendix A. Twist-2 and -3 LCDAs of Light Pseudoscalar Meson 

Two particle twist-2 and twist-3 light-cone distribution amplitudes of light pseudoscalar 



mesons are defined by the following nonlocal matrix elements [ p^ : 

{Pip'MvHqixm = -tfpfip f due'-P'-y+'^^'-'^Mu), (71) 

<^ 



with fp being the decay constant of the light pseudoscalar, fip = Mp/{mi + 1712) (mi and 
1712 are the masses of the constituent quarks in the pseudoscalar), and z = y — x. Here (j){u) 
is the twist-2 light-cone distribution amplitude; (j)p{u) and (paiu) are two-particle twist-3 
distribution amplitudes. The above definitions can be combined into the below nonlocal 
matrix element: 

{P{p')\qa{y)Qpim = ^ f\^-p'-y^^^p'- 

4 JO 

X |/750(t^) - /ip75 (^M^) - W"^^'^) } . (73) 

Neglecting the three-particle twist-3 light-cone wave function, the asymptotic forms of 
the above distribution amplitudes are given as 

(piu) = 6u{l-u), (74) 
Mu) = 1, (75) 
Mu) = 6u{l-u). (76) 

Appendix B. The Evolution of Cj(/i) 
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The renormalization group equation for the Wilson coefficients Cj(/i) is written as follows 



/i|^C(/i) = 7^C(/i) 



(77) 



Here 7 is the anomalous dimension matrix, which can be calculated by the perturbative 
theory and expanded in order of the coupling constants as and aem'- 

7 = irTs + yiTl Ts + — Te + ^^ TsJ H (78) 

471 47r 47r (47r)2 ^ ^ 

The LO anomalous dimension matrix of the above equations has the explicit form 



7, 



(0) 



/ ^ 
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-2 
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-2 


2 
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3Af 
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-6 
N 



















-22 


22 


-4 
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3Af 
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SAT, 










6- ^ 


-6 , 2/ 
TV 3 


-2/ 


2^ 




3Af 
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6 
Af 


-6 








-2/ 


2/ 


-2/ 


-6{JV2-1) 1 


3Af 
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3Ar 


























~-2{u-d/2) 


2{u-d/2) 


-2{u-d/2) 


2(«-d/2) 


3N 
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3N 
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3iV 


3 


3Af 
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V 





~2{u-d/2) 


2{ti-(i/2) 


-2(M-d/2) 


2(«-d/2) 


3N 
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3Af 
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2/ 
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N 
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-6 








6{Af2-l) 
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-6 
N 


6 





6 


N / 



(79) 



where N is the color number, / is the active flavor number, and u and d denote the number 
of the active up- and down-type flavors respectively. 

Appendix C. Some Useful Feynman Parameter Integrals 

In calculation of the perturbative diagrams shown in Fig. 1, one might encounter some 
Feynman parameter integrals which involve nontrivial infrared divergence. To deal with 
the infrared divergence, as mentioned in preceding sections, the dimensional regularization 
(DR) and massive gluon (MG) scheme are applied. Below, we give the explicit calculation 
of some useful Feynman parameter integrals in the above two regularization schemes. 

First, we deal with the integrals in DR scheme. In the DR scheme (here we take d = 4+2a 
and a > 0), the integrals involving the infrared divergence are written as follows: 



1 rl-h I 

dti I dt2-, — TT-j — 

dt'\ I dtoz — ; TT"; 

Jo {tiih + t2u)y-^ 

Jo (-tiM^-" 



1 

u 

1 

u 



+ + -In^M - L12 1 - - 

2 u 



2a 



a 



1 

u 



1 ^ , \nu 
2 + lnMH 

a 1 ~ ""^ 

1 ln(-M)- 2 27 -vr^ 
^ + — + 



(80) 
(81) 



a 



6 



2\n{-u) + -\n^( 



-u 



(82) 
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Here Li2(x) is the dilogarithm function. It is defined by 

^ ln(l -t) 



Li2(a;) = 



t 



-dt . 



(83) 



The Feynman parameter integrals in the MG scheme are hsted as foUows: 



i-ti 



1 



{ti{ti + t2u) + {1 - h - t2)\) 



1 
U 



dti I dt2 
Jo 

1 115 
-In^ X + ln(-u) InA - 2 Iniiln A + - In^ ii - Li2(l ) + -n' 



dti 



i-h 



dU 



to 







1 

u 



dti 



In A - 1 + In M + 



(il(il+i2C^) + (l-il -i2)A) 

In 14 



dig 







1 - u 

{l-ti){l-t2) 
—tit2U + (1 - ti - t2)\ 



(84) 



(85) 



1 
u 



^ In^ A - (In(-ii) - 2) In A - 2 In(-ii) + ^ In^(-ii) + ^ + ^ 



(86) 



When we calculate the above integrals in the MG scheme, the following equations about 
dilogarithm function may be useful: 



(87) 
(88) 
(89) 
(90) 
(91) 



Li2(-a;) + Li2(--) 

X 


6 


-^In^x {x> 0) 


Li2(x) + Li2(-) 

X 


_7r2 
3 


-\n^ X — minx ( 
2 ^ 


Li2(ia;) + Li2( — ) 

X 
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1 1 2 ^1 

— - m X + -nlnx 

Zi Zj 


1 

Li2(-ix) + Li2(-) 

X 
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In^ a; — -tt In x 

2 2 


hi2{x) + Li2(l - x) 
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lnxln(l — x) 



27 



REFERENCES 



[1] M.Wirbel, B,Stech, and M.Bauwe, Z.Phys. C 29, 2637(1985); 34, 103(1987). 

[2] A. Ali and C. Greub, Phys. Rev. D 57, 2996 (1998). 

[3] H.Y. Cheng and B. Tseng, Phys. Rev. D 58, 094005 (1998). 

[4] A. Abada et ah (APE Collaboration), Phys. Lett. B 365, 275 (1996); J.M. Flynn 
and C.T. Sachrajda, Heavy flavours II, 402-452, (|hep-lat/971"0057|) , J.M. Flynn et 
al.(UKQCD Collaboration), Nucl. Phys. B 461, 327 (1996); L. Del Debbio et al. 
(UKQCD Collaboration), Phys. Lett. B 416, 392 (1998). 

[5] P. Ball and V.M. Braun, Phys.Rev.D 55, 5561, (1997); JHEP 9809:005, (1998), (f^- 
ph/ 980239^) . 



[6] A. Khodjamirian, R. Riickl, S. Weinzierl and O. Yakovlev, Phys.Lett.B 410, 275 (1997); 

A. Khodjamirian, R. Riickl and C.W. Winhart, Phys.Rev.D 58, 054013 (1998); A. 

Khodjamirian, R. Ruckl, Heavy flavours II, 345-401, ( |hep-ph / 980 1443|) . 
[7] A. Szczepaniak, E.M. Henley, S.J. Brodsky, Phys. Lett. B 243, 287 (1990). 
[8] Y.Y. Keum, H.N. Li and A.I. Sanda, |hep-ph/000400^ ; Phys.Rev.D 63, 054008 (2001). 
[9] CD. Lii, K. Ukai and M.Z. Yang, Phys.Rev.D 63, 074009 (2001); CD. Lii, M.Z. Yang, 

|hep-ph/ 00112381 . 
[10] CP. Lepage and S.J. Brodsky, Phys. Rev. D 22, 2157 (1980). 

[11] M. Beneke, G. Buchalla, M. Neubert and C.T. Sachrajda, Phys. Rev. Lett. 83, 1914 
(1999). 

[12] H.Y. Cheng, Phys.Lett.B, 335, 428 (1994); Phys.Lett.B, 395, 345 (1997). 
[13] A. Ah, G. Kramer and CD. Lii, Phys.Rev. D 59, 014005 (1999) and Phys.Rev. D 58, 
094009 (1998). 

[14] Y.H. Chen, H.Y. Cheng, B. Tseng and K.C Yang, Phys. Rev. D 60, 094014 (1999). 
[15] J. Chay, Phys.Lett.B 476, 339 (2000). 

[16] M. Beneke, G. Buchalla, M. Neubert and C.T. Sachrajda, Nucl. Phys. B591, 313 (2000). 
[17] D. S. Du, D.S. Yang and CH. Zhu, Phys. Lett. B 488, 46, (2000). 
[18] T. Muta, A. Sugamoto, M.Z. Yang and Y.D. Yang, Phys. Rev. D 62, 094020 (2000). 
[19] M. Beneke, G. Buchalla, M. Neubert and C.T. Sachrajda, |hep-ph/ 0007256 . 



[20] M.Z. Yang and Y.D Yang, Phys. Rev. D 62, 114019 (2000); ^ep-ph/ 00 12208 
[21] X.C He, J.P. Ma, CY. Wu, |hep-ph/0008159| . 
[22] J. Chay, C. Kim, [hep-ph/000921^ . 

[23] H.Y. Cheng and K.C. Yang, Phys. Rev. D 63, 074011, (2001); |hep-ph/0012T5^ . 

[24] D.S. Du, D.S. Yang and CH. Zhu, |hep-ph/00082T^ ; |hep-ph/0l02a77| , to be published 

in Phys. Lett. B. 
[25] M. Beneke, |hep-ph/0009328l submitted to J.Phys.G. 

[26] For a review, see G. Buchalla, A.J. Buras and M.E. Lautenbacher, Rev. Mod. Phys. 

68, 1125 (1996); A.J. Buras, Probing the Standard Model of Particle Interactions, |hep-| 
I ph/980647T| . 

[27] R. Fleischer, Phys. Lett. B 321, 259 (1994); Phys. Lett. B 332, 419 (1994); Z. Phys. C 
62, 81 (1994). 

[28] A.J. Buras and L. Silvestrini, Nucl. Phys. B 548, 293 (1999). 

[29] A.J. Buras, M.Jamin, M.E. Lautenbacher and P.H. Weisz, Nucl. Phys. B 400, 139 
(1993). 

[30] T. Huang, Q. Shen and P. Kroll, Z. Phys. C 50, 139 (1991). 



28 



[31] Y.Y. Keum and H.N. Li, Phys. Rev. D 63, 074006 (2001). 

[32] J.F. Donoghue, E. Golowich, A.A. Petrov and J.M. Scares, Phys. Rev. Lett. 77, 2178 
(1996). 

[33] C.H. Chen and H.N. Li, Phys. Rev. D 63, 014003 (2001). 
[34] A. Khodjamirian, |hep-ph/001227I| . 

[35] V.L. Chernyak and A.R. Zhitnitsky, Phys. Rep. 112, 173 (1983); V.M. Braun and LB. 
Filyanov, Z. Phys. C 48, 239 (1990); P. Ball, J. High Energy Phys. 9901, 010 (1999). 



29 



TABLES 



QCD 
v^oemcients 


H = 5.0 GeV 




fjL = 2.5 GeV 

IVT 


T n 




1.024 + 0.012i 
0.144 - 0.076i 


1.017 
0.188 


1.034 + 0.024i 
0.123 -O.lOOi 


1.037 
0.109 


03 
a| 
al 
as 


0.003 + 0.002i 
-0.027 - O.OMi 

-0.033 - 0.007i 
-0.003 - 0.003i 
-0.036 - 0.012i 
-0.039 - 0.005i 


0.002 
-0.029 

-0.029 
-0.005 
-0.033 
-0.033 


0.004 + 0.004Z 
-0.029 - 0.017i 

-0.036 - 0.007i 
-0.002 - 0.005i 
-0.037 -O.Olli 
-0.040 - 0.004i 


0.004 
-0.040 
-0.040 
-0.010 
-0.040 
-0.040 


07 X 10^ 
r^a^ X 10^ 

r^^Og X 10^ 
ag X 10^ 
a^o X 10^ 
afo X 10^ 


11.9 + 2.8i 
36.8 - lO.Qz 

35.0 - 6.2i 
-936.1 - 13.4i 
-81.8 + 58.8i 
-85.2 + 63.5z 


13.8 
36.8 

36.8 
-928.4 
-141.4 
-141.4 


0.0 + 5.4z 
45.0 - 5.2i 
44.2 + 3.1i 
-953.9 - 24.5i 
-58.3 + 86.1i 
-60.3 + 88.8i 


7.6 
39.8 

39.8 
-957.3 
-74.0 
-74.0 



TABLE I. The QCD coefficients a^(7r7r) at NLO and LO for the renormahzation scales at 
/X = 5 GeV and /x = 2.5 GeV, where = 2m^ / mi,{mu + rnd)- 
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FIGURES 







(a) 



(b) 



(c) 



(d) 







(e) 



(f) 



(g) 



(h) 



FIG. 1. Order of corrections to hard-scattering kernels and T^^ . The upward quark hnes 
represent the ejected quark pairs from b quark weak decays. 





(a) 



(b) 



FIG. 2. Two kinds of topology for penguin contractions. 




Q5(0) 



FIG. 3. An example of the vertex corrections for the operator Q5(0) in coordinate space in the 
case of (hn- insertion. 
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